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ABSTRACT 

We present evidence that the existence of a first order phase transition in com- 
pact U(l) with Wilson action is not related to monopole loops wrapping around 
the toroidal lattice, as has been previously suggested. Our analysis is based on 
the suppression of such loops by 'soft boundary conditions' that correspond to 
an infinitely large chemical potential for the monopoles on the boundary, during 
the updating process. It is observed that the double peak structure characteris- 
tic for the first order phase transition reappears at sufficiently large lattice sizes 
and separations from the lattice boundary. 

1. Introduction 

During this workshop we have heard about plenty of indications that confinement 
in SU(2) and SU(3) gauge theories can be traced back to an underlying U(l) gauge 
structure. The dual Meissner effect mechanism suggests monopoles to condense and 
confine the electric field into a narrow flux tube. It is widely believed — but still 
debated — that the maximal abelian projection uncovers the relevant U(l) degrees of 
freedom in this scenario. 

In 4d compact U(l) lattice gauge theory with Wilson action, monopoles are de- 
fined ab initio as gauge invariant objects, in contrast to non-abelian theories. It 
was demonstrated clearly in ref. [I] that large monopole loops dominate the confining 
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heavy quark potential at large distances in 4d compact U(l) gauge theory. This raises 
the question about the impact of the dynamics of monopole loops onto the critical 
behaviour near the confinement- deconfinement phase transition in U(l). 

Despite the 'simplicity' of the U(l) gauge theory, the order of its phase transition — 
which is relevant to the existence of a continuum limit of the lattice construction — has 
been under debate for 15 years by now, as we can see from table [p. Some years ago, it 

Table 1. Order of the U(l) phase transition. 



author 


max. lattice size 


boundary conditions 


year 


order 


Lautrup et ali 
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Mutter et alJ 
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Jersak et al.u 
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1 st 


Bhanot et al.3 
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Lang et alJI'E 


SHt(lO) 
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Baig et alia 


#4(16) 
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2 nd 


Lippert et al.i 


#4(24) 


wrapping 
loop suppressing 


1994 


1 st 



was argued in ref. |J that the lattice heuristics of metastabilities is greatly affected by 
artifacts from the periodicity on a finite lattice: the wrap-around of large monopole 
loops (around the lattice with toroidal topology) might cause the go-slow of any local 
updating algorithm and thus might fake first order behaviour on a finite size system. 
More recently evidence was presented in support of this view by a study of a lattice 
with trivial homotopy group (without such wrap-arounds! )i@. On a smallish lattice, 
equivalent in volume to a periodic 10 4 lattice, Lang and Neuhaus observed no sign of 
first order behaviour. 

In this contribution we intend to study the anatomy of this system in a variant 
approach. We will be less drastic in changing the lattice topology, remaining on 
a 4d hypercubic lattice. We inhibit wrapping monopole loops during the update 
by suppression of monopole currents to traverse the boundary. This is achieved by 
introduction of an infinitely large chemical potential^ on the boundary, that confines 
monopole currents inside the volume, but allowing sjpin waves to cross. The work 
presented here is an extension of our previous studies^ to a 32 4 lattice. 



2. Monopole suppressing boundary conditions 

Following ref. [n], we start out from the monopole current definition on the lattice, 
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Fig. 1. Plaquettes contributing 
to m M (x). 



Fig. 2. Measurements on nested 
shells (3d projection). 



that is given by the difference between physical flux Q^ u (x) and plaquette flux Q^ u (x). 

As the updating procedure is sequential in terms of the links, we need to clarify 
the interrelation between gauge links and monopole current elements. Let us consider 
the update of a link U p in a 3d cube sitting on the boundary with normalvector e M . 
We search for the plaquettes associated with the Dirac sheets n pa (x), n pu (x) n p - a (x) 
and n p - v {x) that contribute to the monopole current m^(x) into this direction e M . 
These plaquettes are depicted in Fig. [1] with dotted lines. According to Fig. [1] the 
update link U p belongs to four 3d boundary cubes. 

The monopole suppressing boundary condition is injected into the updating pro- 
cedure in the following way: if, for any of these cubes, the change in flux is larger than 
2tt, the proposed link change is rejected in the Monte Carlo. Therefore, the monopole 
current in //-direction in the respective 3d cube is a conserved quantity in Monte 
Carlo time. This allows to enforce configurations with any given number of current 
elements across the boundary. We initialize the simulation with no such crossings at 
all, which implies containment of monopole currents in the volume throughout the 
simulation. 

As the number of monopole currents in a boundary cube is a conserved quantity 
the update rule presented fulfills detailed balance. The procedure of suppressing the 
monopoles can be viewed as being equivalent to an infinitely large chemical potential 
for the monopoles on the cube where the monopole current is suppressed: in the 
equivalent action 5* = Sw + S\, Sw represents the standard Wilson action and 
S\ is the additional term which suppresses monopole currents m M at the lattice slice 
x p = 0, _^ 

S\ = \-Y^ \ m n( x )l A ^ oo. (3) 

Obviously, S\ explicitly breaks translation symmetry. 

In Fig. EL the monople content of a confined configuration taken from a thermalized 
ensemble on an 8 4 lattice at (3 — 1.0 is plotted for illustration. Fig. |3|a shows the 
situation on the standard hypertorus and Fig. |3|b refers to the implementation with 
monopole containment. In Fig. |3|b, no monopole current can traverse the boundaries 
of the lattice. The required absence of traversing currents was monitored throughout 



Fig. 3. a) Monopole cluster wrapping around the standard hyper-toroidal lattice and b) suppression 
of currents across the boundaries leading to containment of the monopole clusters. 



Fig. 4. Plaquette on shells of 16 4 -lattice as 
function of /3.The square region is drawn to 
indicate the size of Fig. @. 



Fig. 5. Plaquette on shells of 24 4 -lattice as 
function of f3. Note the change in scales 
compared to Fig. [|. 



our simulation? 
3. Results 

The simulation was carried out on a sequence of lattices, scanning for the location 
of the phase transition, and in search for first order metastable behaviour. The 
selected lattice sizes were 8 4 , 16 4 , 24 4 , and 32 4 . We performed about 100.000 sweeps 
at each value, that were chosen within the suspected metastable region. We used 
standard Metropolis Monte Carlo updating. 

We measured plaquettes and monopole density, with sampling performed on 
nested shells (see the 3d projection of the shells in Fig. |2|), in order to search gen- 
uine 4d bulk behaviour. Notice that we have traded the unwanted topological effects 
due to wrapping monopole loops against surface effects, i.e. penetration phenomena! 
On a sufficiently large lattice, we would expect with our 'experimental layout' to 
distinguish between first and second order behaviour. 

The results on the 8 4 lattice showed no evidence whatsoever for first order be- 
haviour. On the 16 4 lattice, we scan in j3 through the phase transition. As be seen 
from Fig. f|, the plaquette displays little /3-dependence in the outer shells. For the 
shells with distance greater than 4 from the boundary, however, we discover a rapid 
transient behaviour. This indicates that the boundary layer thickness close to phase 

°In our 4d visualization, see ref. [12], the discrete euclidean time is represented by colors. If an arrow 
traverses the spatial hyperplanes the colors of the arrow's head and tip differ. The boundary in time 
direction is given by the color pink. 



Fig. 6. A (3 scan of the time series on the 8 4 Fig. 7. j3 scan of time series on the outer 24 4 

shell for the 16 4 -lattice. shell for the 24 4 -lattice. 



transition is not prohibitively large for our method to work. Detailed inspection (see 
Fig. H) of the time series on the 8 4 shell with best signal-to-noise ratio, however, does 
not reveal indications for metastabilities. 

We suspect that the 16 4 lattice does not suffice to reveal true bulk behaviour and 
proceed to a 24 4 lattice. The plaquette is plotted across the phase transition in Fig. |5|. 
It shows a much steeper behaviour on the inner shells. Note the change in scale in 
Figs. |] and [5| The region of Fig. |5| transforms to a square drawn in Fig. |j. 

Turning now to the time history, we display the Monte Carlo dynamics as seen 
on the 24 boundary shell (see Fig. |7|), the 16 4 shell (see Fig. |8|), the 10 shell (see 
Fig. 0), and the 4 4 shell (see Fig. |I(]), in a (3 scan across the phase transition. From 
these figures we conclude that (3 = 1.007125 is the most promising setting to study 
the critical dynamics. 

Therefore let us take a closer look at (3 = 1.007125. Fig. O presents the time 
series of the plaquette from the 24 4 down to the 4 4 shell. The fluctuations increase 
as one proceeds to the inner shells, due to a decrease in self averaging. On top of 
those fluctuations, however, the 8 4 shell clearly exhibits long metastabilities of typical 
length of 10000 Monte Carlo sweeps. On the 4 4 shell, the fluctuations become so large 
that the two state metastability signal starts to fade out. 

To further improve the situation we extend our simulation to the 32 4 lattice. 
The resulting pattern of time histories is similar. But the competing short range 
fluctuations are reduced and leave the long range flip-flop behaviour more pronounced. 



8. (3 scan of time series on the 16 4 shell for Fig. 9. (3 scan of time series on the 10 4 shell for 
the 24 4 -lattice. the 24 4 -latticc. 



10. (3 scan of time series on the 4 4 shell for Fig. 11. Scan through different shells of 

the 24 4 -lattice. plaquette time series for the 24 4 -lattice at 

(3 = 1.007125. 



Fig. 12. Histogram of plaquette after smoothing of the series on the 32 4 -lattice at = 1.0090 for 
the 12 4 shell. 



In order to reduce the high frequency modes in the time series, we apply a simple 
filter techniquelil, averaging the series over a sliding window with length I, according 
to P(U) = l/lJ2j=iPi+j-i- The resulting histogram is shown in Fig. [L2|. It pro- 
vides definite evidence of a two-peak structure characteristic for a first order phase 
transition. 

4. References 



1. J. D. Stack and R. J. Wensley, Nucl. Phys. B371 (1992) 597, R. J. Wensley, 
Nucl. Phys. B (Proc. Suppl.) 34 (1994) 204. 

2. B. Lautrup and M. Nauenberg, Phys. Lett. 95B (1980) 63. 

3. K.-H. Mutter and K. Schilling, Nucl. Phys. B200 [FS4] (1982) 362. 

4. J. Jersak, T. Neuhaus and P. M. Zerwas, Phys. Lett. 133B (1983) 103. 

5. Th. Lippert, G. Bhanot, K. Schilling and P. Ueberholz, Nucl. Phys. B (Proc. 
Suppl.) 30 (1993) 912. 

6. C. B. Lang and T. Neuhaus, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 543. 

7. C. B. Lang and T. Neuhaus, Nucl. Phys. B431 (1994) 119. 

8. M. Baig and H. Fort, UAB-FT-338 preprint 1994. 

9. Th. Lippert, A. Bode, V. Bornyakov and K. Schilling, Nucl. Phys. B (Proc. 
Suppl.) 42 (1995) 684. 

10. V. Bornyakov, V. Mitrjushkin and M. Miiller-Preussker, Nucl. Phys. B (Proc. 
Suppl.) 30 (1993) 587. 

11. T. DeGrand and D. Toussaint, Phys. Rev. D22 (1980) 2478. 

12. A. Bode, Th. Lippert, K. Schilling and P. Ueberholz, Int. J. Mod. Phys. C 4 
(1993) 147. 

13. W. Bock, H. G. Evertz, J. Jersak, D. P. Landau, T. Neuhaus and J. L. Xu, 
Phys. Rev. D41 (1990) 2573. 



